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A  previously  reported  [7.  Chem.  Phys.  1991.  95,  106;  Chem.  Phys.  Lett.  1993,  206,  137]  Fourier  transform 
method  for  computation  of  classical  intramolecular  mode-to-mode  energy-transfer  rate  coefficients  is  extended 
to  four-atom  mdecules.  HONO  and  C2H2  are  used  as  test  cases.  The  method  involves  the  Fourier  transform 
of  the  time  variation  of  a  local-mode  bond  energy  for  an  ensemble  of  trajectories.  A  two-mode,  coUinear  model 
is  employed  to  demonstrate  that  the  transform  is  expected  to  contain  a  series  of  spectral  bands  at  frequencies 
corresponding  to  the  mode-to-mode  energy-transfer  rates.  Heavy-atim  blocking  and  constrained  motion  methods 
are  employed  to  determine  the  individual  band  assignments.  The  results  for  both  HONO  and  C2H2  are  in  good 
accord  with  the  total  relaxation  rate  extracted  from  decay  plots  of  the  local-mode  energy. 
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I.  Introductioa 

We  have  previously  described'’^  a  Fourier  transform  (FT) 
method  for  extracting  classical  intramolecular  vibrational  energy- 
transfer  rates  (FTIVR).  The  FTIVR  method  involves  the 
computation  of  the  time  dependence  of  some  local-mode  bond 
energy,  £b(0,  in  the  molecule  from  the  ensemble  average  of  a  set 
of  classical  trajectories  computed  using  standard  techniques.^ 
This  bond  energy  exhibits  characteristic  fluctuations  due  to 
energy-transfer  processes  between  the  various  modes  in  the 
molecule.  Consequently,  the  FT  of  £(,(/)  should,  in  principle, 
yield  a  spectrum  that  contains  bands  at  frequencies  corresponding 
to  the  mode-to-mode  rate  coefficients. 

In  previous  applications,  we  have  investigated  mode-to-mode 
energy-transfer  rates  in  the  ethynyl  (C2H)  radical'  and  disilane.^ 
In  the  case  of  C2H,  three  distinct  bands  appear  in  the  FTIVR 
spectrumin  the  frequency  range 0.0 1-0. 18  ps~'.  The  frequences 
of  these  bands  are  interpreted  to  represent  the  average  mode- 
to-mode  energy-transfer  rates  for  IVR  via  the  stretch/stretch 
and  two  stretch/bend  pathways.  It  is  shown  that  these  mode- 
to-mode  rates  are  in  good  accord  with  the  total  relaxation  rate 
extraaed  from  decay  plots  of  the  C~H  local-mode  energy.  We 
have  also  applied  the  FTIVR  method  to  intramolecular  eneigy 
transfer  in  ditilane.2  Theresultisaseriesofbroadpeaksspanning 
the  range  0.01-0.3  ps~'  with  a  global  maximum  at  0.07  ps~'.  This 
suggests  that  the  IVR  rates  in  disilane  are  described  by  a  set  of 
coefficients  whose  values  lie  in  the  range  0.01-0.3  ps*'  with  a 
most  probable  value  of  0.07  ps*' .  This  interpretation  was  shown 
to  be  consistent  with  the  total  relaxation  rate  of  the  Si-H  local 
mode  extracted  by  least-squares  fitting  of  the  decay  curves 
obtained  from  the  trajectory  calculations.  Unfortunately,  the 
resolution  of  the  spectral  bands  was  insufficient  to  permit  the 
individual  coeHicients  to  be  determined.  This  difficulty  arises, 
in  part,  because  of  the  large  number  (153)  of  mode-to-mode 
coefficients  involved. 

In  the  present  paper,  we  first  explore  the  FTTVR  method  in 
a  simple  model  system  which  permitsa  more  detailed  investigation 
than  is  possible  in  more  complex  molecular  systems.  We  then 
show  that  the  FTIVR  method  can  be  effectively  extended  to  four- 
atom  molecules.  HONO  and  acetylene  are  chosen  as  test  cases. 
This  choice  is  dictated  by  the  availability  of  accurate  intramo¬ 
lecular  potentials  and  by  the  fact  that  C2H2  and  HONO  represent 
both  linear  and  nonlinear  cases. 

n.  AppBcatioa  to  a  Smpie  Model  System 

Consider  an  ensemble  of  N  local  viluational  modes.  If  the 
coupling  terms  between  these  modes  are  zero,  each  will  correspond 
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FigaK  1.  Definition  of  variables  in  the  two-mode,  coUinear  model. 

to  a  good  action  variable,  and  the  corresponding  mode  energies 
will  be  constants  of  the  motion.  A  Fourier  transform  of  the  mode 
energies  would  therefore  exhibit  only  a  single  component  at  zero 
frequency.  If,  however,  the  coujding  terms  are  nonzero,  in¬ 
tramolecular  energy  flow  between  the  modes  will  be  observed 
and  the  mode  energies  will  vary  with  time.  The  basic  pronise 
of  the  FTIVR  method  is  that  siiKX  the  fluctuations  of  the  mode 
energies  are  produced  by  mode-to-mode  energy  transfer,  the 
corresponding  rates  for  these  underlying  processes  must  be 
embedded  in  the  time-varying  behavior  of  the  local-mode  energies. 

Consider  a  simple  coUinear  system  consisting  of  two  barmmiic 
oscUlators  with  a  quadratic  coupling  potential.  A  diagram  of 
thissystemisshowninFigurel.  WetakethesystemHamilumian 
to  be 

H  =  +  F2V2/W2  +  + 

[Pi^/Tm^  +  /»4V2«4  +  *2(^2  -  ■Re2)'l  +  ^('•3  -  ( J) 

where  the  m,  (<  =  1,  2,  3,  4)  are  the  atomic  masses.  Pi  are  the 
momenta,  and  the  are  interatomic  distances  asdefined  in  Figure 
1.  The  local  mode  energies  are 


£,  =  [F,V2m,  +  +  A:,(r,  -  (2) 

and 

=  [FjV2m3  +  P//2in4  +  kj(rj  -  R^)^]  (3) 

If  the  coui^ing  comtant  A  is  zero,  £|  and  £2  wiU  be  constants  of 
the  motion.  For  A  0,  the  modes  wiU  be  coupled  and  energy 
transfer  will  occur  between  them. 

We  have  investigated  the  nature  of  the  energy  transfer  for  the 
above  system  for  a  total  energy  of  O.SO  eV  initiaUy  contained  in 
mode  1.  The  system  parameters  are  given  in  Table  1.  Figure 
2,  a  and  b,  shows  the  time  variation  of  £|  for  A  equal  to  0.45  and 
0.60  eV/A^  respectively.  Several  qualitative  features  of  the 
energy  transfer  between  inodes  1  and  2  are  immediately  evident. 
First,  the  mode-to-mode  energy  transfer  produces  characteristic 
oscillations  in  the  local-mode  energy.  Second,  the  osdllatioiis 
resulting  from  energy  transfer  between  modes  1  and  2  have  higher 
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Flgere  2.  Energy  in  local  mode  1  as  a  function  of  time:  (A)  X  =  0.4S 
eV/A^(B)X  =  0.60eV/A^.  Time  is  given  in  units  tu,  where  1  tu=  1.019 
X  10-'<  s. 


FiguK  3.  Fourier  transform  of  the  local  mode  bond  energy  £i(()  given 
in  Figure  2A. 


TABLE  1:  ParaoKteis  of  the  Two-Mode,  CoOinear  IVR 
Model 


ki 

1.00  eV/A* 

1.00  A 

«2 

10.00  amu 

1.00  eV/ A* 

1.50  A 

mj 

10.00  amu 

1.00  A 

mi 

10.00  amu 

nu 

10.00  amu 

frequencyoscillationssuperimposeduponthem.  Third, the mode- 
to-mode  energy  transfer  rate  increases  as  the  coupling  constant 
X  increases. 

Figure  3  shows  the  Fourier  transform  of  £|  given  in  Figure  2a 
for  X  =  0.45  eV /  A^.  The  FT  spectrum  shows  three  features.  The 
most  important  of  these  is  a  band  at  65.S  cm~’,  which  corresponds 
to  a  frequency  of  2.0  X  10**  s*'  and  a  period  of  5.0  x  10"**  s.  In 
the  motoailar  time  units  being  used  in  Figure  2a.b.  this  corresponds 
to  a  period  of  49.07  tu  (1  tu  =  1.019  X  ICh**  s).  An  inspection 
of  Figure  2a  shows  that  after  500  tu  the  local-mode  energy  £| 
hasexecutedabout9.5completeoscil]ations.  Obviously, the maj<xr 
FT  spectral  band  at  65.5  enr*  corresponds  to  the  frequency  at 
which  energy  is  transferred  between  the  modes.  A  similar  Fourier 
transform  of  Fgure  2b  yields  a  three-band  spectrum  with  the 
major  peak  at  85.2  cm~',  which  corresponds  to  the  mode-to-mode 
energy-transfer  rate  seen  with  X  =  0.60  eV/A*. 

The  small  peak  in  the  FT  spectrum  at  176.9  cm~‘  is  the  result 
of  the  higher  frequency  component  seen  in  Fgure  2a.  These 
osdllations  reflect  energy  transfer  in  and  out  of  the  coupling 
term  in  eq  1  rather  than  actual  transfer  to  nude  2.  This  is  easily 
seenfromthetransfonnofthecouidingtermalone.  Thisisshown 
in  Figure  4  for  the  case  X  ^  0.45  eV/A*.  The  only  prmninent 
feature  is  a  band  at  176.9  cm~'.  Finally,  the  remaining  band  in 


1/cm 

rigmc4.  Fourier  transfer  oftbe  time  variatianofthe  coupling  term  with 
X  «  0.45  eV/A*. 


Rgare  5.  Dependence  of  the  frequency  of  the  energy-transfer  FT  band 
on  the  coupling  constant,  X. 


0.02  0.04  0.06  0.08  0.10  0.12 
1/<Half-Ufe>  (1/t.u.) 

Rgnted.  Dependence  of  the  frequency  of  the  energy-transfer  FT  band 
on  the  inverse  of  the  average  half-life,  computed  for  different 
values  of  X. 

the  FT  spectrum  of  £|(r)  is  a  peak  at  1 1 1 .4  cm~' .  This  frequency 
is  the  difference  between  the  energy-transfer  band  at  65.5  cm~* 
and  the  coupling  term  band  at  176.9  cm'*.  The  peak  at  111.4 
cm~'  may  therefore  be  identified  as  the  comlnnation  band  (vpc- 
vsi)- 

Since  the  frequency  of  the  energy-transfer  band  in  Fgure  3 
measures  the  rate  of  mode-to-mode  energy  transfer,  we  would 
expect  it  correlate  directly  with  the  magnitude  of  the  coupling 
term  connecting  the  two  modes.  Figure  5  shows  the  calculated 
frequency  of  the  energy-transfer  band  as  a  function  of  the  coupling 
constant  X.  The  direct  correlation  is  obvious. 

The  data  shown  in  Figure  2a,b  make  it  clear  that  energy  transfer 
in  the  model  system  cannot  be  described  by  a  first-order  rate  law. 
It  is  therefore  not  possible  to  deflne  a  rate  coeffldent  in  the  usual 
manner  as  the  negative  slope  of  a  logarithmic  decay  curve.  We 
can,  however,  obtain  a  second  measure  of  the  tracer  rate  by 
computation  of  an  average  half-life,  <ri/2>.  This  is  most 
conveniently  done  by  averaging  the  results  of  M  trajectories 
differing  only  in  the  initial  phase  of  local  mode  1  and  computing 
the  times  required  for  the  peak  energies  in  mode  1  to  deoease 
to  half  of  th^  initial  values,  (ti/j)  is  defined  as  the  ensemble 
average  of  these  values.  Since  <ri/2}~'isameasureoftbeeneigy- 
transfer  rate,  we  wmild  expect  a  direct  correlation  to  exist  between 
the  values  of  <li/2>''  obtained  for  diflerent  values  of  X  and  the 
frequences  of  the  energy-transfer  band  in  the  FT  spectrum  of 
£i.  figure  6  shows  that  the  cmrelation  between  these  quantities 
is  excellent. 
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Figac7.  Heavy-atom  Mocking  effects.  (A)  Time  variation  of  the  energy 
in  local  mode  2  for  £  =  O.SO  eV  initially  paititioaed  into  mode  1  with 
Mi  -  1000  amu  and  X  =  0.43  eW/k^.  Other  parameters  are  as  given 
in  Table!.  (B)  Fourier  transform  of  local  mode  energy  l.£i(r),forthe 
case  described  in  Figure  7A. 

If  the  mode  1  mode  2  path  is  effectively  blocked  by  assigning 
a  very  large  mass  to  atom  3,  we  expect  the  energy-transfer  rate 
to  approach  zero  and  for  the  FT  band  corresponding  to  this  mode- 
to-mode  rate  to  decrease  significantly  in  intensity  or  to  vanish 
entirely.  Figure  7a  shows  Ezit)  for  the  case  M3  ^  1000  amu  and 
X = 0.45  eV /  A^,  with  the  remaining  parameters  being  those  given 
in  Table  1.  Obviously,  the  energy-transfer  path  is  essentially 
blocked.  No  more  than  10%  of  the  toul  energy  initially  present 
m  mode  1  ever  finds  its  way  into  mode  2.  Figure  7b  shows  the 
FT  spectrum  of  £|({)  for  this  case.  As  can  be  seen,  the  energy- 
transfer  band  at  6S.Scm~'  has  vanished  alongwiththeoombination 
band  at  111.4  cm''.  Only  the  mode  coupling  band  remains.  This 
peak  is  slightly  red-shifti  to  163.8  cm*'  due  to  the  effect  of  the 
heavy-atom  mass. 

The  above  analysis  shows  that  the  frequency  of  the  energy- 
transfer  band  in  the  FT  spectrum  correlates  directly  with  ti 
magnitude  of  the  mode-to-mode  coupling  term  and  with  other 
measures  of  the  transfer  rate  such  as  <(i/2)~'.  Examination  of 
the  time  dependence  of  £|  shown  in  Figure  2a  demonstrates  that 
the  FT  band  frequency  corresponds  exactly  to  the  energy-transfer 
rate.  Elimination  of  energy  transfer  by  heavy-atom  blocking 
eliminates  the  energy-transfer  band  entirely.  It  is  therefore 
appropriate  to  take  the  frequency  of  the  FT  energy-transfer  band 
to  be  the  characteristic  mode-to-mode  transfer  rate. 

We  now  consider  how  the  existence  of  many  vibrational  modes 
in  a  three-dimensional  molecular  system  would  be  expected  to 
alter  the  results  obtained  for  the  two-mode,  coilinear  model  shown 
inFigurel.  First,  there  will  be  many  additional  energy-transfer 
pathways,  each  with  its  own  characteristic  frequency  and  coupling 
terms.  We  will  therefore  expect  the  FT  spectrum  to  contain 
many  bands.  Some  of  these  will  correspond  to  mode-to-mode 
energy-transfer  rates;  others  will  be  assodated  with  rates  of  energy 
transfer  to  and  from  the  coupling  terms.  In  addition,  we  expect 
to  find  combination  bands  and  perhaps  some  overtone  bands. 
Identification  of  these  bands  will  therefore  pose  a  significant 
problem  which  is  analogous  to  the  experimental  problem  of 
assigning  IR  absorption  bands  to  particular  vibrational  modes  in 
the  molecule. 

Because  of  the  presence  of  many  transfer  pathways,  we  expect 
entropy  effects  to  significantly  reduce  the  ampUtude  of  the 
oscillations  in  the  local-mode  energies.  However,  we  still  expect 


them  to  be  present.  If  there  are  enough  coupled  modes,  the 
osciUation  amplitudes  may  be  reduced  to  the  point  that  the 
envelope  function  begins  to  resemble  an  exponential  decay  which 
will  permit  us  to  extract  a  total  energy-transfer  coefTicient  out 
of  mode  1,  K'r,  which  will  be  the  sum  of  all  of  the  mode-to-mode 
rates, 

J 

In  such  a  case,  we  expect  the  sum  of  the  frequencies  of  the  FT 
energy-transfer  bands  to  correlate  directly  with  ATV  determined 
fromtheslopeofthelogarithmicdecaycurve.  We  have  previously 
shown  that  this  is  the  case  for  C2H  and  disilane.  In  the  present 
paper,  we  also  show  that  it  is  true  for  HONO. 

Finally,  since  there  will  be  many  and  varied  coupling  terms 
present  in  the  complex  potential  representing  a  real  molecule 
with  numerous  vibrational  modes,  we  expect  the  time  variation 
of  the  local-mode  energy  to  exhibit  many  high-frequency 
oscillations  superimposed  upon  lower  frequency  oscillations  which 
correspond  to  the  energy-transfer  bands  whose  amplitudes  will 
be  a  decreasing  function  of  time. 

The  calculation  of  an  FTIVR  spectrum  has  the  dual  advantages 
that  it  is  very  easy  to  compute  and  it  contains  all  of  the  mode- 
to-mode  rates.  Tlie  two  major  diRiculties  associated  with  the 
method  involve  band  resolution  and  assignment.  If  there  are 
many  mode-to-mode  bands,  it  may  be  impossible  to  resolve  them. 
This  is  analogous  to  the  problem  of  resolving  the  rotational  bands 
in  a  large  molecule.  The  first  difficulty  may  be  minimized  to 
some  extent  by  extending  the  computation  time  for  the  trajectories 
which  will  increase  the  FT  resolution.  The  second  problem  may 
be  attacked  using  heavy-atom  blocking  and  constrained  motion 
methods.  These  methods  are  described  in  more  detail  in  the 
following  section. 


A.  HONO.  Thepotentialenergysurfaceusedinthecalculation 
is  that  employed  in  our  previous  energy-transfer  studies.*  Bond 
stretching  potentials  are  represented  by  a  summation  of  Morse 
potentials. 


I 

where  the  summation  runs  over  all  bending  motions  of  the 
molecule.  The  torsional  potential  is  given  by 

Fi  =  ]^a,cos(iV)  (7) 

I 

where  the  £,  are  bond  lengths,  61  are  bending  angles,  and  9  is  the 
dihedral  an^e. 

B.  Acetylet.  The  analytical  potential  energy  surface  devel¬ 
oped  by  White  and  Schatz^  is  employed  for  acetylene.  Thissurface 
is  based  on  fits  to  ab  initio  calcuhtions  of  the  molecular  force 
field,  experimentally  derived  energies  of  formation,  and  the 
isomerization  barrier.  The  fitting  procedure  used  for  C2H 
combines  a  LEPS  potential  with  a  three-body  Sorbie-Murrell 
function*  to  fit  the  C2H  force  field,  the  isomerization  barrier,  and 
the  C2H  dissociation  energies.  For  C2H2,  two  C2H  fragment 
potentials  are  combined  with  an  empirical  methylene  potmitial 
and  a  four-body  Sorbie-Murrell  function*  to  fit  the  acetylene 
force  fkid,  the  vinylidene  minimum  energy,andothCT  information. 


m.  Potential  Energy  Sarfiaces  and  Nanwrkal  Procednres 


Quadratic  functions  are  used  for  the  bending  motions 
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TABLE  2:  C2H2  and  HONO  Noraial-Mode  Frequencies 

molecule  mode  no. 

mode  description 

freq  (cm"') 

CjH,*  1 

C-H  stFctch 

3400 

2 

C-H  stretch 

3314 

3 

C-C  stretch 

2022 

4-5 

HCCH  bend* 

797 

6-7 

HCCH  bend* 

569 

HONO  1 

O-H  stretch 

3438 

2 

N— 0  stretch 

1640 

3 

O-N  stretch 

908 

4 

HON  bend 

1306 

5 

ONObend 

540 

7 

torsion 

649 

*  Computed  on  the  White-Schatz  potential  surface,  ref  S.  *  Doubly 
degenerate  vibrational  mode. 


The  fundamenul  vibrational  frequencies  for  ail  of  the  HONO 
and  acetylene  normal  modes  computed  using  the  above  potentials 
are  summarized  in  Table  2. 

C.  Niunerical  Procedures.  The  local  bond  mode  energy  is 
deflned  as 

£b(0  =  (^(0V2Mb)  +  nK(0]  (8) 

where  the  reduced  mass  fib,  momentum  PtO),  and  bond  length 
rb(()  are  those  appropriate  to  the  bond  in  question.  The  bond 
potential,  Kb[rb(f)]>  taken  to  be  a  Morse  function  with 
parameters  chosen  so  that  it  corresponds  closely  to  the  bond 
potential  on  the  full  global  surface. 

Hamilton’s  equations  of  motion  are  integrated  using  a  fourth- 
order  Runge-Kutta  integrator.  The  acetylene  intramolecular 
dynamics  are  computed  for  a  time  period  196.6  ps  using  a  fixed 
integration  step  size  of  5.0  X  10-5  ps.  Conservation  of  energy  to 
six  significant  figures  is  generally  obtained.  For  HONO,  the 
internal  dynamics  are  followed  for  82  ps  using  a  step  size  of  l.S 
X  10~*  ps.  In  this  case,  energy  conservation  was  generally  three 
significant  digits  or  better. 

Initial-state  selection  for  acetylene  is  accomplished  by  inserting 
the  zero-point  vibrational  energy  (ZPE)  computed  from  the 
potential  field  into  each  of  the  normal  modes  with  the  molecules 
in  their  equilibrium  configurations.  Subsequently,  the  equations 
of  motions  are  integrated  for  a  random  peri^  of  time  after  which 
a  selected  amount  of  additional  energy  is  inserted  into  the  local 
C-H  stretching  mode  in  the  form  of  kinetic  energy.  For  HONO, 
the  ZPE  of  each  normal  mode  is  partitioned  between  kinetic  and 
potential  energy  with  appropriate  averaging  over  the  phases  of 
the  normal  coordinates,  followed  by  selective  excitation  of  the 
O-H  stretch  local  mode.  The  coordinates  and  momenta  are  then 
scaled  to  the  desired  total  energy.  The  details  of  this  procedure 
have  been  fully  described  elsewhere.^ 

Ensemble  spectra  are  computed  by  averaging  the  Fourier 
transforms  of  each  of  40  trajectories  that  differ  only  in  the  random 
period  or  phase  averaging  used  in  the  initial-state  selection. 

Molecules  containing  four  atoms,  such  as  HONO  and  C2H], 
havesix  vibrational  modes  and  1 5  mode-to-mode  rate  coefficients. 
We  therefore  expect,  in  principle,  a  FTIVR  spectrum  containing 
as  many  as  1 S  energy-transfer  peaks  plus  higher  frequency  bands 
resulting  from  energy  transfer  to  and  from  the  coupling  terms 
along  with  overtone  and  combination  bands.  The  assignment  of 
thesespectral  features  to  particular  mode-to-roode  energy-transfer 
processes  is  generally  a  diHlcult  problem.  We  have  found  two 
techniques,  heavy-atom  blocking  and  constrained  motion,  to  be 
particular  useful  methods  to  unravel  the  FTIVR  spectrum. 

Heavy-atom  blocking  is  a  technique  originally  developed  by 
Rabinovitch  et  al.  in  their  energy  flow  experiments  on  cyclic, 
organometallic  compounds.*  In  these  experiments,  it  was  found 
that  if  two  carbon  rings  are  separated  by  a  heavy  metal  atom, 
energy  flow  between  the  rings  is  severely  restricted.  The  results 
obtained  with  the  two-mode,  collinear  model  serve  as  a  simple 


Time  (ps) 


FigB«  8.  (a)  Ensemble  average  decay  curve  for  the  O-H  local-mode 
energy  in  HONO.  The  plot  is  the  ensemble  average  of  40  trajectories. 
The  initial  states  have  ZPE  in  all  of  the  vibrational  modes  and  sufEdent 
exdtation  of  the  O-H  local  mode  to  yield  a  total  energy  of  2.31  eV.  (b) 
FTIVR  spectrum  obtained  from  the  average  of  the  Fourier  transforms 
of  the  local  O-H  bond  energy  for  each  of  the  trajectories  used  to  obtain 
the  decay  plot  shown  in  Figure  la.  The  flgure  inset  showing  HONO  with 
all  atoms  written  in  nonboidface  type  indicates  normal  masses  arc  used 
for  all  atoms. 

example  of  the  effect  The  same  principle  may  be  employed 
effectively  to  obtain  the  band  assignments  in  the  FTIVR  spectrum. 
For  example,  in  the  acetylene  molecule,  H~|C— jC—H,  if  the 
normal  mass  of  jC,  12  amu,  is  replaced  by  a  heavy  mass  of  100 
amu,  the  energy  pathway  from  the  H~|C  bond  to  the  2C-H 
stretch  motion  is  effectively  blocked.  As  a  result,  the  spectral 
bands  associated  with  the  blocked  energy-transfer  pathway(s) 
will  be  either  missing  or  significantly  reduced  in  intensity.  A 
comparison  of  the  FTIVR  spectra  obtained  with  and  without 
heavy-atom  blocking  will  therefme  permit  the  missing  or  intensity- 
reduced  peaks  to  be  assigned  to  the  blocked  pathways. 

The  method  of  constrained  motion  involves  the  omission  of  ail 
initial  energy  in  certain  modes  of  the  molecule.  This  omission 
serves  to  signiflcantiy  reduce  or  eliminate  the  coupling  between 
those  modes  and  others  in  the  molecule.  As  a  result,  the  spectral 
bands  associated  with  energy  transfer  between  those  modes  will 
be  reduced  in  intensity  or  missing.  This  intensity  change  permits 
the  bands  assignments  to  be  made.  For  example,  if  the  initial 
zero-point  torsional  mode  energy  is  omitted  in  HONO,  there  will 
be  no  initial  momentum  components  popendicular  to  the 
mdecular  plane.  In  principle,  the  system  will  be  constrained  to 
planar  motion.  Consequently,  the  importance  of  the  torsion 
pathway  in  I VR  will  be  substantially  reduced,  and  those  spectral 
bands  associated  with  energy-transfer  pathways  involv^  the 
torsional  mode  will  decrease  in  intensity. 

IV.  Results  and  Discnnoo 

Figure  8a  shows  the  decay  of  log(£b(0]  fot  HONO  at  a  total 
energy  2.31  eV,  corresponding  to  zero-point  energy  in  all  normal 
modes  with  the  excess  energy  initially  present  as  local-mode 
exciution  in  the  O-H  bond.  The  decay  curve  is  obtained  by 
averaging  over  the  results  of  40  trajectories.  The  general  form 
of  the  decay  curve  is  as  anticipated.  There  is  an  a|^>roximateiy 
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linear  decay  with  several  low-frequency  oadlladons  superimposed. 
Additional  higher  frequency  components  are  also  present. 

A  total  relaxation  rate  coefficient  out  of  the  C-H  mode,  K-r, 
can  be  extracted  from  these  data  by  fitting  the  decay  curve  to  a 
first-order  rate  law: 

£b(»)  =  £b(»“0)«P(-^Tf)  (9) 

so  that 

log(£b(01  »  log[£b('=0)l  -  ( 10) 

By  obtaining  a  linear  least-squares  fit  of  eq  10  to  the  first  portion 
of  the  data  shown  in  Figure  8a,  we  obtain  ATt  1.86  X  10"  s~'. 
Since  we  expect  K^  to  be  the  sum  of  all  of  the  mode-to-mode  rate 
coefficients  out  of  the  local  O-H  mode,  it  is  clear  that  the  individual 
mode-to-mode  rate  coefficients  should  lie  in  the  range  0.01-0. 1 86 
ps'',where0.01  ps~'  is  an  assumed  lower  limit  for  such  coefficients. 

The  ensemble  averages  of  the  Fourier  transforms  of  each  of 
the  40  trajectories  yield  the  results  shown  in  Figure  8b.  The 
energy  present  in  each  mode  varies  as  vibrational  relaxation 
proce^.  It  is  therefore  expected  that  each  mode-to-mode 
coefficient  will  appear  as  a  broad  band  covering  a  range  of 
frequencies  that  is  a  reflection  of  the  energy  range  sampled  in  the 
trajectories.  This  can  be  seen  in  the  spectra  presented  in  Figures 
8-10.  Several  broad  bands  are  present  in  the  spectral  region 
between  1  and  4  cm*'.  The  summation  of  these  bands  results  in 
six  distinct  peaks. 

The  frequencies  at  which  these  bands  occur  represent  the 
average  m^e-to-mode  rate  coefficients.  The  only  remaining 
problem  is  to  determine  the  modes  associated  with  each  peak 
appearing  in  Figure  8b.  The  band  assignments  are  made  using 
heavy-atom  blocking  and  constrained  motion  methods. 

Figure  9a  shows  the  FTIVR  spectrum  obtained  if  the  masses 
of  the  two  oxygen  atoms  and  nitrogen  in  HONO  are  assigned  a 
value  of  100  amu.  This  assignment  will  effectively  block  the 
energy  flow  from  the  O-H  stretch  to  all  other  modes  except  the 
HON  bend  and  out-of-plane  torsional  mode.  The  ensemble 
FTIVR  spectrum  obtain^  with  this  mass  combination,  using  the 
same  random  seed,  exhibits  three  peaks,  marked  as  1. 2,  and  3, 
which  appear  at  the  same  frequencies  as  the  correspondingly 
marked  bands  in  Figure  8b.  Bands  4,  S,  and  7  are  missing. 
Consequently,  they  must  correspond  to  energy-transfer  pathways 
blocked  by  the  iinvy  atoms. 

Intramolecular  vibrational  energy  redistribution  in  HONO 
has  been  studied  extensively  by  Thompson  and  co-workers.^  They 
have  concluded  that  the  major  energy-transfer  pathway  from  the 
O-H  stretch  is  via  direct  transfer  to  the  HON  bend  followed  by 
transfer  from  the  HON  bend  to  the  torsional  mode.  We  therefore 
assign  peak  1,  the  lowest  frequency  spectral  band  in  Figures  8b 
and  9a,  to  the  O-H  stretch/torsion.  Peak  2,  which  corresponds 
to  the  second  fastest  energy-transfer  rate,  is  assigned  to  the  HON 
bend/t(xsk>n,  and  peak  3,  the  fastest  rate,  is  assigned  to  the  O-H 
stretch/HON  bend. 

Figure  9b  shows  the  ensemble  FIT  VR  spectrum  obtained  when 
the  nitrogen  and  oxygen  atoms  in  the  N-0  end  group  are  both 
assigned  masses  of  100  amu.  We  would  expect  this  mass 
combination  to  block  energy-transfer  pathways  involving  the 
terminal  N-O  stretch.  As  expected,  bands  1,  2,  and  3  are  still 
presmit  in  the  spectrum.  In  addition,  three  peaks  denoted  4, 5. 
andOarealsoevident.  Thesenewpeaksareexpectedtoberelated 
to  the  energy  flow  paths  Mocked  by  the  heavy  oxygen  atom  in 
the  H-O-N  moiety  in  Figure  9a.  That  is,  we  expect  peaks  4,  S, 
and  6  to  involve  the  O-N  stretch  and  the  O-N-O  bend. 

By  employing  the  constrained  motion  method,  we  may  obtain 
the  additional  information  required  to  make  the  band  assignments. 
The  FTIVR  spectrum  given  in  Figure  9c  is  that  obtained  when 
the  initial  ZPE  for  the  ONO  bend  is  <»iitted  and  the  terminal 
N-O  atoms  are  both  assigned  masses  of  100  amu.  Comparison 


FigurcP.  (a)  FTIVR  spectrum  of  HONO  with  ZPE  plus  sofliciem  O-H 
local-mode  excitatkm  to  yMd  a  total  internal  energy  of  2.31  eV.  The 
atoms  in  the  0-44-0  mtnety  are  each  assigned  a  heavy-atimi  mass  of  1 00 
amu.  This  is  denoted  by  the  boldface  type  in  the  figure  inset.  (b)FTTVR 
spectrum  oi  HONO  under  the  same  conditions  as  in  Figure  2a  except 
only  the  terminal  N  and  O  atoms  have  a  heavy  mass  of  100  amu.  Othm 
arenonnal.  (c)  FTIVR  spectrum  of  HONO  under  the  same  conditions 
as  Figure  2a  except  the  initial  ZPE  of  the  ONO  bending  mode  is  set  to 
zero,  (d)  FTIVR  spectrum  of  HONO  under  same  conditions  as  tbosc 
given  for  Figure  8  b  except  the  ZPE  of  the  torsional  mode  is  initially  set 
to  zero.  All  masses  are  normal  as  indicated  by  the  figure  inset 
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Figore  10.  (a)  FTIVR  spectra  for  acetylene  obtained  from  the  ensemble 
average  of  the  Fourier  transforms  of  each  of  40  trajectories.  The  initial 
states  have  ZPE  in  all  of  the  vibrational  modes  and  suflicient  excitation 
of  the  C-H  local  mode  to  yield  a  total  energy  of  3.18  eV.  Dotted 
spectrum:  normal  mass  for  all  atoms.  Uppersolid-linespecttum:  terminal 
(nonexcited)  C  and  H  atoms  have  mass  100  amu.  Others  have  normal 
mass.  Lower  solid-line  spectrum:  atoms  in  the  C-C-H  moiety  have 
mass  100  amu.  (b)  Ensemble  average  (40  trajectories)  FTTVR  spectra 
foracetylene.  Thedottedspectrumisthesameasuppersolid-linespectrum 
shownin(a).  Thesolid-linespectrumisobtainedundercooditionsidentical 
to  those  for  the  dotted-line  spectrum  except  the  initial  ZPE  of  the  CCH 
bending  mode  is  set  to  zero. 

of  parts  b  and  c  of  Figure  9  shows  that  peak  S  has  disappeared 
while  peaks  4  and  6  remain.  We  therefore  assign  band  S  to  the 
O-H  stretch/ONO  bend  energy-transfer  pathway. 

The  remaining  assignment  is  made  by  noting  that  peak  6  does 
not  appear  in  the  FTIVR  spectrum  sho^  in  Figure  8b  where  all 
frequencies  corresponding  to  the  rates  of  major  energy-transfer 
patte  should  be  present.  We  also  note  that  the  frequency  of  peak 
6  is  exactly  double  that  of  peak  4.  We  therefore  take  peak  6  to 
be  an  overtone  of  peak  4  and  assign  the  O-H  stretch /O-N  stretch 
rate  to  the  frequency  of  peak  4. 

A  consistency  cross-check  may  be  obtained  by  setting  the  initial 
ZPE  of  the  HONO  torsional  mode  to  zero  with  all  atoms  assigned 
their  actual  mass.  Under  these  conditions,  there  are  no  initial 
momentum  components  perpendicular  to  the  molecular  plane. 
Therefore,  in  principle,  the  system  will  be  constrained  to  planar 
motion,  and  energy-transfer  pathways  involving  the  bnsknal  mode 
will  be  blocked.  [We  note  that  because  of  numerical  integration 
inaccuracies,  some  torskmal  motion  may  actually  occur  in  the 
calculations.  However,  its  magnitude  sho^d  be  greatly  reduced.] 
Figure  9d  shows  the  resv  dng  FflYR  ensemUe  spectrum. 
Comparison  with  Figure  8b  shows  that  peak  2  is  eithn  absent 


TABLE  3:  Mode-to-Mode  Rate  Cocfficteats  for  Interaal 


Energy  Trunfcr  in  CiHi  and  HONO 

molecule 

peak  no. 

transfer  path 

K{pr') 

C,H, 

1 

C-H  stietcb/CCH  bead 

0.026 

2 

C-H  stretch/ torsion 

0.041 

3 

HCC  bend/totrion 

0.0s  1 

4 

C-H  stretch/C-C  stretch 

0.061 

5 

C-H  stretch/HCC  bend 

0.071 

6 

C-H  stretch/C-H  stretch 

0.082 

HONO 

I 

O-H  stretcb/torsioD 

0.0060 

2 

HON  bend/torsion 

0.024 

3 

O-H  stretcb/HON  bend 

0.049 

4 

O-H  stretch/O-N  stretch 

0.037 

5 

O-H  stretch/ONO  bend 

0.061 

7 

O-N  stretch/N-O  stretch 

0.08S 

or  significantly  reduced  in  magnitude.  This  is  consistent  with 
our  previous  assignment  of  band  2  to  the  HON  bend/torsicm 
pathway.  Finally,  we  note  that  peak  7  only  appears  in  Figures 
8b  and  where  all  atoms  have  their  normal  mass.  This  suggests 
that  this  band  is  associated  with  the  O-N  stretch/N-O  stretch 
path.  All  of  the  mode-to-mode  rate  coefficients  obtained  from 
the  FTIVR  spectra  for  HONO  are  summarized  in  Table  3. 

A  second  cross-check  on  the  accuracy  of  the  FTIVR  results 
may  be  obtained  by  using  them  to  compute  the  total  relaxation 
rate  of  the  O-H  niode  and  comparing  the  value  obtained  with 
the  IVR  relaxation  rate  computed  from  eq  10  and  the  data  in 
Figure  8.  Since  relaxation  of  the  O-H  mode  can,  in  principle, 
occur  by  energy  transfer  to  the  O-N  stretch,  the  HON  bend,  the 
ONO  bend,  the  torsion  mode,  and  N-O  stretch,  we  would  expect 
the  total  relaxation  rate  coefficient,  K^,  to  be  given  by 

+  (11) 

where  Kt  is  the  frequency  corresponding  to  peak  i.  The  HONO 
potential  energy  surface  used  in  the  present  study  contains  no 
direct  coupling  terms  between  the  O-H  and  N-^  stretching 
modes.  Consequently, weexpectAN-otobezeroforthispotential. 
The  fact  that  we  have  not  bwn  able  to  identify  an  FTIVR  band 
corresponding  to  the  O-H  stretch/N-O  stretch  pathway  is 
therefore  not  surprising.  Since  Kj  is  obtained  directly  from  the 
trajectory  results,  while  the  Kt  are  obtained  from  the  FTIVR 
band  frequencies,  the  extent  to  which  eq  1 1  holds  constitutes  an 
internal  check  mi  the  consistency  of  the  results  obtained  by  the 
two  methods. 

From  Table  3  for  HONO,  we  obtain  (JSTi  +  JiTs  +  A4  +  JSC5)  = 
1.S3  X  10"  s~'.  The  spectral  resolution  for  each  band  is  0.20 
cm"'  or  0.06  X  10"  s-'.  We  therefore  have  (Ai  -F  ATj  +  A4  +  Ki) 
=  (1.53  ±  0.24)  X  10"  r'.  The  least-squares  fitting  of  eq  10  to 
the  data  shown  in  Figure  8  yields  Ky  —  1.86  x  10"  s"'.  The 
major  uncertainty  in  this  result  is  related  to  the  arbitrary  choice 
of  the  initial  time  period  over  which  the  least-squares  fitting  is 
performed.  By  varying  this  choice  within  reasonable  limits,  we 
estimate  that  Kj  may  be  varied  by  about  ±10%.  Thus,  the 
consistency  check  compares  (1.S3  ±  0.24)  X  10"  s~'  with  (1.86 
±  0.19)  X  10"  s~'.  Since  the  error  limits  of  these  values  have 
significant  overlap,  we  conclude  that,  within  the  limitations 
imposed  by  the  spectral  resolution  and  the  uncertainty  present 
in  the  analysis  of  the  trajectory  results,  the  mode-to-i^e  rate 
coefficients  obtained  from  the  FTIVR  spectra  and  the  trajectory- 
computed  overall  O-H  relaxation  rate  are  consistent. 

As  an  example  of  a  linear  molecule,  we  have  api^ed  the  FTIVR 
method  to  obtain  an  analysis  of  intramolecular  energy  flow  in 
acetylene.  The  spectra  shown  in  Hgure  10a  cmrespo^  to  the 
ensemble  average  of  the  Fourier  transform  of  each  of  40 
trajectmies  for  acetylene  at  a  total  energy  of  3.18  eV.  This  energy 
comprises  zero-point  energy  in  all  normal  modes  with  the  excess 
energy  initially  present  as  local-mode  excitation  of  one  of  the 
H-C  bonds.  The  dotted  line  shows  the  FTIVR  spectrum  for 
normal  CjHz  while  the  solid  lines  present  results  for  cases  involving 
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heavy-atom  blocking.  The  top  line  is  jbtained  when  the  terminal 
C  and  H  atoms  in  the  C-H  bond  not  initially  excited  are  assigned 
a  mass  of  100  amu.  Comparison  with  the  dotted  line  shows  that 
peak  6  is  missing  while  t^  rest  of  the  bands  remain.  Peak  6  is 
therefore  assigned  to  the  H-C  stretch/C-H  stretch  pathway. 

The  bottom  solid  curve  in  Figure  10a  shows  the  result  when 
all  atoms  in  the  C-C-H  moiety  are  assigned  a  mass  of  100  amu. 
Under  these  conditions,  we  expect  only  three  energy-transfer 
pathways  to  be  open,  the  H-C  stretch/HCC  bend,  the  H-C 
stretch/torsion,  and  the  HCC  bend/torsion.  Comparison  with 
the  dotted  curve  in  Figure  10a  shows  that  peaks  1 , 4,  and  6  have 
vanished.  Only  peab  2,  3,  and  S  remain.  On  the  basis  of 
arguments  similar  to  those  used  for  HONO  and  in  the  analysis 
of  the  FTIVR  spectrum  for  the  ethynyl  radical,'  we  assign  the 
highest  frequency  band  S  to  the  H-C  stretch/HCC  bend  and 
peak  3  to  the  HCC  bend/torsion.  The  lower  frequency  band  2 
is  interpeted  to  result  from  the  C-H  stretch/torsion  pathway. 

By  elimination  and  by  examination  of  the  predicted  rates,  band 
4  and  1  are  assigned  to  the  C-H  stretch/C-C  stretch  and  H-C 
stretch /CCH  bend,  respectively.  The  numerous  bands  at  higher 
frequencies  either  are  associated  with  energy  transfer  in  and  out 
of  the  numerous  coupling  terms  in  the  potential  as  the  model 
calculations  suggest  or  correspond  to  overtone  or  combination 
bands. 

Figure  10b  shows  the  FTIVR  spectra  for  CiHj  with  two  heavy 
atoms,  C  and  H.  The  dotted  curve  is  the  same  as  the  top  solid 
curve  in  Figure  10a.  The  solid  curve  in  Figure  10b  is  the  ITIVR 
spectrum  obtained  with  the  ZPE  for  the  CCH  bend  and  torsional 
modes  initially  set  to  zero.  It  is  seen  that  bands  1 , 2,  and  3,  which 
are  associated  with  the  H-C  stretch/  "^CH  bend,  the  H-C  stretch/ 
torsion,  and  the  HCC  bend /torsion,  respectively,  are  considerably 
reduced  as  we  would  expect  from  the  principle  of  constrained 
motion.  This  suggests  that  the  band  assignments  are  correct. 
The  average  individual  mode-to-mode  rate  coefTicients  obtained 
from  the  FTIVR  spectra  for  C2H2  are  summarized  in  Table  3. 

V.  Swnmary  and  CoacliBimis 

By  employing  the  Fourier  transform  of  a  local-mode  bond 
energy,  mode-to-mode  rate  coefficients  for  classical  intramolecular 
vibrational  relaxation  have  been  obtained  for  HONO  and  C2H2. 
For  both  HONO  and  C2H2,  six  distinct  peaks  appear  in  the  FTIVR 


spectrum.  The  frequencies  at  which  these  bands  appear  are 
interpreted  to  correspond  to  the  average  mode-to-mode  rate 
coefficients  for  IVR.  To  determine  the  modes  associated  with 
each  individual  peak,  heavy-atom  bioclong  and  constrained  motion 
methods  have  been  used  to  close  certain  energy-transfer  paths 
during  vibrational  mode  relaxation.  Comparison  of  the  resulting 
FTIVR  spectra  with  those  obtained  under  normal  conditions 
permits  the  accurate  assignment  of  individual  bands. 

The  present  results  combined  with  earlier  studies  on  the  ethynyl 
radical'  show  that  this  FTIVR  method  works  well  on  three-  and 
four-atom  molecules.  Previous  studies^  involving  an  eight-atom 
molecule,  disilane,  demonstrate  that  the  method  permits  the 
average  value  of  the  mode-to-mode  rate  coefTicients  to  be 
extracted,  but  limitation  on  spectral  resolution  has,  to  date, 
prevented  the  determination  of  individual  coefficients  in  such 
large  systems. 
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